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ABSTRACT 

We examine the prospects of making a joint analysis of neutrino oscillation at two 
baselines with neutrino superbeams. Assuming narrow band superbeams and a 100 kt 
water Cerenkov calorimeter, we calculate the event rates and sensitivities to the matter 
effect, the signs of the neutrino mass differences, the CP phase and the mixing angle 
^^13. Taking into account all possible experimental errors under general consideration, 
we explored the optimum cases of narrow band beam to measure the matter effect and 
the CP violation effect at all baselines up to 3000 km. We then focus on two specific 
baselines, a long baseline of 300 km and a very long baseline of 2100 km, and analyze 
their joint capabilities. We found that the joint analysis can offer extra leverage to 
resolve some of the ambiguities that are associated with the measurement at a single 
baseline. 



1 Introduction 

Although the existing data from the Super-Kamiokande experiment and various other 
corroborating experiments offer very strong indications of neutrino oscillations, the ap- 
pearance experiment, i.e., the appearance of a flavor different from the original one, has 
not been convincingly performed. If neutrinos indeed oscillate, the oscillation parame- 
ters, including the leptonic CP phase, have to be determined with sufficient accuracy. 
Furthermore, the well-known MSW matter effect has to be tested by experiments. 
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In spite of the various ongoing and planed neutrino oscillation experiments, additional 
experiments with very long baseline are needed, at least for the test of the matter ef- 
fect. The recently approved superbeam facility 0], which will be available towards the 
later part of this decade, offers the possibility of a very long baseline (VLBL) experi- 
ment which, in conjunction with other oscillation experiments, can test thoroughly all 
properties of neutrino oscillations. 

Among all neutrino oscillation experiments, the long baseline (LBL) experiments are 
particularly attractive. Since the neutrino beams are produced in an accelerator accord- 
ing to definite physics criteria with the detector site chosen accordingly, the experiment 
can be conducted in a more controlled fashion to maximize the physics output. Hence 
the LBL experiments will allow us to make detailed analyses of the oscillation param- 
eters so as to provide a complete picture of the physics of neutrino oscillation. As one 
example of such experiments, a project called H2B is under discussion 0, |^, |^. The 
neutrino super-beam for H2B would be from the newly approved high intensity 50 GeV 
proton synchrotron in Japan called HIPA and the detector, tentatively called the 
Beijing Astrophysics and Neutrino Detector (BAND), is envisioned to be a 100 kt water 
Cerenkov calorimeter (WCC) with resistive plate chambers (RPC) [0 located in Beijing, 
China. The distance from HIPA to Beijing is about 2100 km. Such a very long baseline 
experiment would be complementary to the recently proposed J2K experiment P] which 
will also use the neutrino beam from HIPA but with the Super-Kamiokande detector or 
its update. The distance from HIPA to Super-Kamiokande is about 300 km. 

In this article, we will examine the prospects of investigating neutrino oscillations at 
H2B in conjunction with J2K so that the joint data at the two widely different baselines 
can be used in a complementary way to provide strong leverage to eliminate some of 
the ambiguities in the determination of oscillation parameters. The joint analysis can 
expand the capability of the parameter search that are not attainable by either of the 
experiments alone. The two baselines can work at their respective favorable energy 
ranges. The present work is to demonstrate this possibility. But we have not search for 
the best narrow beam energies for the two baselines. Assuming a narrow band meson 
beam and the above mentioned 100 kt WCC with RPC, we simulate the event rates for 5- 
year operation. The sensitivity of the event rates for the various oscillation parameters 
will be explored. The present work can be regarded partly as a continuation of the 
study of H2B Refs. 0, @, |^ and an initial exploration of the idea of joint analyses of 
two detectors which we think is appropriate for oscillation physics. In Sec. 0, we discuss 
some of the fundamentals of neutrino oscillation and LBL experiments. In Sec. ^ we 
present some of our numerical results. We present the joint analyses of the data of two 
detectors in Sec. 4. Finally, in Sec. 5, we present our conclusions. 
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2 Fundamentals of neutrino oscillation and LBL ex- 
periments 

If we accept all current data, there will be three distinctive mass scales provided by the 
five categories of experiments: long baseline, short baseline accelerator experiments such 
as LSND, atmospheric, solar, and reactor. If the LSND data are excluded, the three SM 
neutrino fiavors are sufficient and no extension of the number of neutrinos beyond that 
of the standard model is necessary. In view of the uncertainty of the LSND data, our 
discussion will be restricted to the 3-fiavor scenario. 

The oscillation of the 3-fiavor neutrinos is a system with a limited number of degrees 
of freedom. The system consists of 2 mass square differences (MSD), three mixing 
angles and one measurable CP phase. These parameters together with the matter effect 
determine the various survival and appearance probabilities 0. The unitary mixing 
matrix in vacuum is generally parameterized as 

(C12C13 C13S12 §1^ 

— C23S12 — C12S13S23 C12C23 — S12S13S23 C13S23 
S12S23 - C12C23S13 -C12S23 - C23S12S13 C13C23 

where sjk = sm{9jk), cjk = cos{9jk), and Sjk = sm{9jk)e'^^ , 9jk defined for j < k is 
the mixing angle of mass eigenstates Uj and Uk, and 6 is the CP phase angle. The 
three mass eigenvalues are denoted as mi, m2, and m^. The two independent MSD are 
Am.21 = m2 — m^ and Amgj = mg — m^. 

In LBL experiments the neutrino beam has to go through matter which gives rise to 
the well-known MSW effect 0. A widely used model for the Earth, called the prelimi- 
nary reference Earth model PREM, is given in [0] and the earth density profile can be 
found in ||rT|. Since for a VLBL experiment the matter density can vary significantly 



along the path of the neutrino beam, in our calculation we perform numerical integration 
of the Schrodinger equation for a realistic treatment of the distance dependent matter 
density. 

The detection of a given neutrino fiavor is through its accompanying charged lepton 
produced by the charge current interaction of the neutrino with the nucleons in the 
detector mass. For a neutrino energy E^y, which is small compared to the mass of 
the W and Z bosons but large enough so that quasi-elastic effect is small, the charge 
current cross sections are given by a^^ = 0.67 x 10~^^cm^i?,^(GeV) for electron and muon 
neutrinos, and apN = 0.34 x 10~^^cm^ii^;y(GeV) for electron and muon anti-neutrinos. 
For the tau neutrino, the above expression is subject to a threshold suppression. The 
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threshold for the production of the tau is Et = rrij- + = 3.46 GeV. A fit of u-j- to 
cross section as a function of the neutrino energy in terms of the ratio of two quadratic 
polynomials can be found in Ref. The signal events of fiavor f3, i.e., the number of 
charged lepton of fiavor f3, from a neutrino beam of fiavor a, to be observed at a baseline 
L is given by 

N,= ^E,,L)a{E,)Po,^p{E,,L)dE,, (2) 
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where ^{Ei,, L) is the total neutrino flux spectrum including the detector size and run- 
ning time period, Pa^/^^E^, L) is the oscillation probability, a{E^) the neutrino charge 
current cross section, and -^Max and -Emin are the maximum and the minimum energies 
of the beam. 

In a narrow band beam the neutrino flux is distributed below a given energy -Epeak- 
The intensity is peaked at -Epcak and decreases rapidly below £'peak- The wide band 
beam contains neutrinos with energy spread out in a signiflcant range of energy. In our 
calculation we will use the realistic beam energies and proflles provided in |T^. Some 
of the narrow band beams together with the wide beam are plotted in Fig. 1. Here 
dNcc/dEy = ^{E^, L)a{E^) is the energy distribution of the charged-current events Ncc 
for one year operation of a 100 kt detector at L=2100 km. 

Since in oscillation experiments, especially in the case of electron neutrino appear- 
ance, the statistics are generally not large. Therefore the error is an important factor in 
the physics extraction. W use the approach of Ref. to estimate the possible statistical 
and systematic errors and to gain a sense of the goodness of the flt. For the electron 
counting experiments the errors and uncertainties arise from the following sources: 

(i) The statistical error in the measurement of the charge lepton of flavor f] which is 
as usual ^/Ns + N^. Nb is the number of measured background events and can be 
expressed as 



(ii) The systematic uncertainty in the calculation of the number of background events, 
which can be denoted as V/sNi,. 

(iii) The systematic uncertainty in the beam flux and the cross section which we denote 



The total error is the quadrature of all these uncertainties. In our calculation we will 
take rp = 0.1, gp = 0.05, and fp = 0.01. 

3 Numerical results for individual baselines 

Presently there are sizable errors in all the oscillation parameters. However, we envisage 
that at the H2B time, Amgg, Am|^, ^23? and 612 will be fairly accurately determined. 
So we will not assign any speciflc errors to them. We focus our investigations on the 
following parameters and effects: matter, MSD sign, CP violation, and ^13. 

3.1 Inputs 

We present numerical results of a 5-year operation with a water Cerenkov detector. The 
detector size is assumed to be 100 kt for all baselines. Sizes other than 100 kt will be 
labeled whenever used. 




(3) 



as gpN, 
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The inputs of the mixing angles and MSB's are from solar, atmospheric and CHOOZ 
experiments. For definiteness we take sin^(26'i2) = 0.8 and sin^(26'23) = 1.0. In most 
of our results we use sin^(2^^i3) = 0.05 for illustration and effects of larger and smaller 
values of 6113, 0.01 < sin2(26'i3) < 0.1, will be investigated. The inputs of MSB Am^^ 
and Am32 are respectively given by Am^^j = 5 x 10~^ eV^ and Am^^.^^ = 3 x 10"^ eV^. 

Presently the sign of the MSB's are unknown so there are 4 possibilities: 



II III IV 

Am 



^32 + + - ^ ^ 



Ami 



After showing the effects of all four sign combinations in the electron event numbers we 
will choose the sign I for illustration. 

3.2 Matter effects 

In Tables 1 and 2 we show the —>■ v^. event rates with and without matter effects for 
a narrow band beam with -Epeak = 4 GeV for both baselines. It is clear that for both 
narrow band and wide band beams the matter effect is significant on electron event 
number at L=2100 km, but negligible at L=300 km. As expected, the z/^ and z/^ events 
show very little matter effect at either distance. The event rates at both baselines can 
be increased if different narrow band beams are used. For example, for L=2100 km the 
-E'peak=6 GeV beam has twice as many electron events as the -Epeak = 4 GeV beam. 

In order to look for the optimum beam energy to measure matter effects at a given 
baseline, we have examined the following ratio, which is approximately the statistical 
significance of the matter effect and is referred to in Ref. as the figure of merit, 

A^e I with matter A^e | without matter / _ \ 

-^matter — a"1\7 " \ J 



AN, 



e 



Here AA"e is the total error of the electron event number, as discussed at the end of Sec. 
2, without the matter effect. Figure 2 shows -Rmatter versus the baseline up to 3000 km 
for several narrow band beams for the four MSB signs combinations. We see that for 
L=2100 km the optimal narrow band beams for the matter effect are with peak energies 
in the range of 4 ~ 6 GeV. For example, as shown in Fig. 2 for the MSB sign I, the 
optimal narrow band beam has the peak energy around -Epeak = 4 GeV. For L=300 km, 
as expected, there is very little statistical sensitivity to the matter effect at all available 
energies. 

Given a narrow band beam with -Epeak = 4 GeV for L=2100 km and -Epeak = 0.7 GeV 
for L=300 km. Fig. 3 shows the electron event rate versus the CP phase with or without 
matter effect. We see that for ^13 to have a fixed value or small range of uncertainties 
the matter effect is experimentally measurable for L=2100 km but hardly observable for 
L=300 km. However in the currently fully allowed range of 613, sin^(20i3) < 0.1, it is 
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even difficult for the 2100 km baseline to distinguish the matter effect from the vacuum 
for the following fact: Since the electron event rate is proportional to sin^(2^i3), the 
electron event rates for sin^(26'i3) = 0.03 with matter effect and for sin^(26'i3) = 0.1 in 
the case of vacuum are the same as can be inferred from Fig. 3, it is not possible to 
distinguish the two. This ambiguity will be reinforced when the error is not negligible. 

3.3 MSD sign effects 

The sensitivity of the event rate to the sign of MSD for sin^(26'i3) — 0.05 is also shown in 
Tables 1 and 2 for E^eak — 4 GeV and 5 = for both baselines, and in Fig. 4 for different 
energies for the two baselines as functions of the CP phase. Tables 1 and 2 show that 
the electron event rates are sensitive to the sign of MSD at the 2100 km baseline. It is 
also interesting to note that for L=300 km there is sensitivity in distinguishing signs I 
and IV in which both MSD are positive or negative from signs II and III in which one 
is positive and the other negative. This general feature is vahd for other values of ^13 
once it is determined. 

In Fig. 4, in which we take sin^(2^i3) = 0.05, it shows clearly that for L=2100 km 
I and II are well separated from III and IV for all values of CP phase. Hence the 
sign of Am|2 should be readily determined with moderate amount of electron neutrino 
appearance data. However, the separation of I from II depends on the value of the 
CP phase. In the region of small, intermediate and large value of the CP phase, the 
sign of Am^Qi can be determined, but around 6 = 130° and 6 = 280° I and II are not 
distinguishable. The signs HI and IV are almost inseparable in the whole region of S. 
Hence the sign of Am^-^ will be very hard to determine if < 0. Then the anti- 

neutrino beam is needed for the determination. For L=300 km. Fig. 4 shows that it is 
difficult to distinguish I, II, HI and IV except in very special values of the CP phase. 

Unfortunately, the above result is only true if ^13 is already known. Similar to the 
situation discussed at the end of the preceding subsection, the significant uncertainty in 
sin^(26'i3) muddies the water. As sin^ (26*13) decreases the electron event rate will also 
be reduced. Therefore, it is difficult to distinguish the signs I and II of small ^13 with 
signs HI and IV with a larger ^13. We demonstrate the decrease of the Icpton event 
rate with sin^(2^i3) in Fig. 4. Hence when the full range of current uncertainty of 613 is 
include, i.e., sin^(26'i3) < 0.1, the sensitivity in distinguishing the MSD sign is lost for 
both basehnes. 

3.4 CP violation effects 

Figures 3 and 4 show the electron event number versus the CP phase, modulo the 
matter effect. The typical total errors are also shown. The dominant error is found to 
be statistical, i.e., from the source (i) as described at the end of Sec. 2. We see that 
although the event rate varies significantly with the CP phase, as the electron event rate 
is not a single valued function of the CP phase, it is ambiguous to determine 6 from the 
electron event number even for a fixed value of 6*13. The caveat of the uncertainty in 6*13 
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discussed in the two previous subsections made the ambiguity even more serious. 

The sensitivity of the electron event rate to the CP phase depends on the beam 
energy as shown in Fig. 5. At some of the beam energies, e.g., 2 and 10 GeV for L=2100 
km and 0.7 GeV for L=300 km, the curves are quite flat, indicating a poor sensitivity 
to the CP phase at such beam energies. Furthermore at almost no energies that one can 
determine a unique CP phase from the electron event number at either 300 km or 2100 
km. 

To investigate the sensitivity we define two ratios involving the two CP conserving 
phases: 5 = 0° and 5 = 180°: 

= AiV,(180°) ' 

where Ne{6), iVe(0°), and A^e(180°) are respectively the electron event numbers for CP 
phases 5, 0° and 180°, and ANe{0°) and AA^e(180°) are the total error at 5 = 0° and 
6 = 180°. We can now defined the figure of merit 0, i.e., the goodness of the fit, for the 
CP violation measurement as the smaller in magnitude of the two ratios: 



CP 



In Fig. 6 we plot Fcp(5) versus the peak energy of the narrow band beam, separately 
for L=2100 and 300 km. We show six values of 5=0°, 30°, 60°, 90°, 120°, and 150°. The 
curves satisfy approximately the relation Fcp{l80° + 6) ^ —Fcp{S). Hence the curves 
for 5 =180°, 210°, 240°, 270°, 300°, and 330° can be inferred as the negatives of the 
above corresponding curves of 5 less than 180°. The left panel is for the 100 kt detector 
and the right panel shows the results for a 1000 kt detector. We see that for the 100 kt 
detector at both baselines the effects of the finite CP phases are within la from each 
other, including the CP conserving case. If we increase the detector size to 1000 kt, the 
CP violation effects can reach to the 2a level for the beams around E^eak — 3-4 GeV 
and 6-7 GeV for 6 = 60°-120° and 240°-300° at L=2100 km, and around Ep^ak ^0.7 
GeV for the similar 6 ranges at L=300 km. 



3.5 Effects of the uncertainty of sin^(2^i3) 

In all the above results we have used sin^(26'i3) = 0.05. Since z/g is proportional to 
sin^(2^i3), the latter is a sensitive parameter for the electron event number. Accordingly, 
the counting experiment of the electron event number may provide a good measurement 
for the value of sin^(2^i3). 

In Fig. 7 we present the electron event number versus the CP phase for different 
sin^(26'i3) values. The error bars indicate the size of the estimated total errors. From 
the total errors, we see how precisely the sin^(2^i3) value can be measured. For example, 
for L=2100 km the curve of sin^(26'i3) = 0.08(0.06) lies about 1.5cr (30") away from that 
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of sin^(2^i3) = 0.1. Then it is difficult to distinguish 0.1 from 0.08 all along the curves. 
Furthermore without knowing the CP phase, it may be difficult to distinguishing 0.1 
at one CP phase to 0.6 at another CP phase. This ambiguity is even more serious for 
L=300 km because there is more variation of the event number as a function of the CP 
phase. 

4 Joint analysis of baselines 2100 and 300 km 

We imagine that major efforts of the very long baseline experiments such as H2B are the 
conffimation of the matter effect, the determination of the MSD signs, the CP phase, 
and ^13. However, there exist difficulties in finding unique solutions for them, given 
the measured electron event rate, as demonstrated in the preceding section. We have 
discussed repeatedly in the previous section the ambiguities caused by the current wide 
range of uncertainty in ^^13. There are other ambiguities which are caused by the multi- 
valueness of the oscillation probability as a function of the oscillation parameters and 
the possibility of overlapping parameter regions. To illustrate the latter ambiguity let 
us consider Fig. 4. For the simplicity of argument, let us ignore any possible errors. 
Suppose a measurement of the electron event rate is 60 at 300 km baseline for a narrow 
band beam with peak energy 0.7 GeV. Then the CP phase can be either around 0° or 
150° for sin^(2^^i3) = 0.05. Similarly, suppose a measurement at the 2100 km baseline 
gives, say, the electron event rate is 40 at 4 GeV. Then CP phase can be either 150° or 
300° for sin^(2^i3) = 0.05. Further, since the value of sin^(2^i3) is unkown, we in fact 
obtain a curve in the 5 — sin^(2^i3) plane for a given electron event number, as shown 
in Fig. 8. Hence the measurement from only one experiment, either at L=300 km or at 
L=2100 km, is not enough to determine CP phase or the value of sin^(2^i3). 

To illustrate the advantage of the joint analysis of two widely different baselines, we 
plot in Fig. 8 sin^(2^i3) vs 6 for measured electron event rates for both 300 km and 
2100 km basehnes at respectively 60 and 40 events for the MSD sign I. In the absence 
of any errors, the intersect of the curves gives unique values of both sin^(2^^i3) and S. 
In reality the situation will be more complicated due to the presence of errors of the 
measurements, and hence the intersect of the two curves will cover a sizable area of the 
sin^(2^i3) vs 5 plane. However, this example shows the possibility of extra leverages one 
can gain with two different baselines. 

In this section we present some of our analyses of such joint measurements, taking 
the advantage of superbeams like HIPA, which can offer multiple narrow band beams of 
different energies. We use different energies at the two baselines. We will plot 2100 km 
baseline vs 300 baseline by simultaneously looking at two different parameters. 

4.1 sin2(2^i3) and the CP phase 6 

In Fig. 9 we show electron event number at L=2100 km versus those at L=300 km for 
fix MSD sign I. Each curve has a fixed value of sin^(2^i3) with the CP phase 6 varies 
in the full possible range from 0° to 360°. The 6 — 0° point is marked by a sohd dot 
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and the S = 180° point by a cross. The direction of increasing S is indicated by the 
arrow on the curve. The curves are generally ellipses and the eccentricity of the ellipse 
is determined by the specific beam energies of the two baselines. 

We fix 0.7 GeV for the 300 km baseline and allow the energy at 2100 km to change. 
The upper diagram of Fig. 9 is at 4 GeV for 2100 km. When sin^(2^i3) increases the 
ellipse moves towards the upper right, i.e., increasing the electron event rate for both 
baselines. This is expected from the fact that the oscillation probability z/g is 

proportional to sin^ (26*13). Since the ellipses of neighboring values of sin^(26'i3) overlap 
significantly, the value of 5 and sin^( 26*13 ) can not be determined uniquely, reflecting 
again the ambiguities discussed in the preceding section. However there are energies at 
which the overlap of the ellipses is minimized. The lower diagram of Fig. 9 shows that 
the ellipses of constant ^^13 are collapsed into lines when the beam energy of the 2100 
km baseline is 6.3 GeV. So in principle the joint measurement allow us to narrow down 
the allowed range of sin^ (26*13). For the fines each measurement still allows two values 
of 5. But the two values of S which fall on top of one another on the line segment will 
be separated when the line becomes an ellipse. So measurements at both 6.3 and 4 GeV 
will offer a better possibility to determined the values of sin^(2^i3) and 6 simultaneously. 

In Table 3 we present, for the case of MSD sign I, some Ep^ak values in GeV of 
narrow band beams where the ellipses of Ne(300) versus Ne(2100) as the CP Phase 
varies from 0° to 360° collapse into fines. At these energies the curves for MSD sign II 
arc ellipses of high eccentricities which approximate lines. For MSD signs III and IV, 
and in the absence of matter cfi^cct the curves are ellipse of very high eccentricities. For 
these energies the combined measurements of electron event at L=2100 km and L=300 
km can provide better measurement for the sin^(2^i3). 

4.2 MSD sign and the CP phase S 

In Fig. 10 we present similar results, but for different MSD signs with fixed sin^ (26*13) = 
0.05. The results without the matter effect are also plotted, with the dotted curves 
denoting MSD sign II or III and dashed ones I or IV. In the absence of the matter effect 
MSD signs I and IV give the same results, so do the MSD signs II and III, as already 
shown in Tables 1 and 2. For the almost overlapped curves of MSD sign III and IV with 
matter effects, the sofid ones denote III and dotted ones IV. 

It is clear from Fig. 10 that in the lower diagram, i.e., 6.3 GeV for the 2100 km 
basefine, it is quite easy to differentiate MSD signs I and II from III and IV, and from 
the case without the matter effect. To make better measurements it is again better to 
take measurements with the line together with the ellipse. 

5 Conclusion 

In the above study of the event rates and the sensitivity to various oscillation parameters 
investigated, we found: 
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(a) At the distance L=2100 km, a narrow band beam with peak energy of about 6 
GeV is optimum for measuring CP violation effects and about 5 GeV for measuring 
matter effects. 

(b) To measure the CP violation effect at a shorter distance such as L=300 km, a 
narrow band beam with lower peak energy (~ 0.7 GeV) is preferable. But the 
matter effect is hardly observable at such a shorter baseline. 

(c) The two baselines, 300 km and 2100 km, are complementary to each other. Through 
the joint analysis of the two baselines, some of the ambiguities associated with the 
measurement at either baselines may be resolved. 

With the optimum narrow band beam, a 5-year operation of a 100 kt water Cerenkov 
detector at a very long distance such as L=2100 km has the following physics prospects: 

(1) The matter effects can be observed. 

(2) The sign of Am|2 niay be determined. 

(3) The sign of Am2]^ may be determined only in favorable situations. 

(4) Evidence exceeding 2-0" of a CP violating phase may be seen in favorable cases for 
a detector size of 1000 kt or with a much longer running time. 

(5) Combined with the analyses of L=300 km, the parameter sin^(2^i3) may be mea- 
sured and the matter effects are more clearly determined. 

In this article we have focused on the — > z/g exclusively. The investigation of the r 
appearance and the inclusion of the beam option in the analysis, which is needed in 
the cases of MSD signs III and IV, i.e., Am32 < 0, will be taken for a future investigation. 
There we will also make a more complete search for the best energies of the two baselines 
for the various parameters. 

We finally note that the statistics are generally low in all the cases discussed. Running 
with higher energy narrow band beam will increase the statistics. However, that may 
be disfavored by the figure of merit (signal to error ratio). Another way to increase the 
statistics is to increase the detector mass. It has been pointed out, however, that there 
is a saturation problem |jl3| caused by the systematic errors which are of the form of 
the errors of types (ii) and (iii) as discussed at the end of Sec. 2. These errors increase 
linearly as the number of events rather than the square root of the number of events as 
is in the case of the statistical error. Hence, when the mass of the detector is increased 
so that the number of events becomes sufficiently large, the systematical error becomes 
dominant. After that, further increase of the detector size may no longer be beneficial. 
In Fig. 11 we show the ratio of AA'^e to Ng as a function of the detector mass. We see 
that according to our general error estimate the best AN^ to ratio can be attained 
is 6%. When the detector reaches 1000 kt the benefit of further increasing the detector 
size is no long significant. 
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Table 1: Event rates of 5-year operation with (without) matter effects for different 
MSD sign choices for a narrow band beam of -Epeak = 4 GeV. The CP-phase is taken to 
be zero. 





electron ^ 


muon ^ 


tau # 


I 


34 (10) 


430 (435) 


10 (11) 


L=2100 km II 


46 (16) 


405 (415) 


11 (11) 


III 


3(16) 


413 (415) 


12 (11) 


IV 


3(10) 


427 (435) 


11 (11) 


I 


159 (157) 


39408 (39407) 


72 (72) 


L=300 km II 


119 (116) 


39535 (39535) 


71 (71) 


III 


114 (116) 


39535 (39535) 


71 (71) 


IV 


154 (157) 


39408 (39407) 


72 (72) 



Table 2: Same as Table 1, but for a wide band beam. 





electron 7^ 


muon 7^ 


tau # 


I 


151 (96) 


2313 (2311) 


448 (453) 


L=2100 km II 


151 (90) 


2326 (2333) 


443 (449) 


III 


39 (90) 


2335 (2333) 


454 (449) 


IV 


49 (96) 


2308 (2311) 


458 (453) 


I 


453 (443) 


271536 (271535) 


731 (731) 


L=300 km II 


359 (348) 


271842 (271842) 


718 (718) 


III 


337 (348) 


271843 (271842) 


718 (718) 


IV 


431 (443) 


271535 (271535) 


731 (731) 



Table 3: Some -Epeak values (GeV) of narrow band beams where the ellipses of Ne(300) 
versus Ne(2100) as CP Phase varies from 0° to 360° collapse into line segments. The 
MSD sign is assumed to be case I. 



^peak(300) 


^peak(2100) 


0.70 


0.750 1.215 1.85 2.30 6.30 


0.80 


0.820 1.10 1.98 2.25 7.60 


0.85 


0.820 1.20 2.05 2.19 8.30 
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Figure 1: The energy E,^ distribution of charged-current events Ncc for one year oper- 
ation of a 100 kt detector. 
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Figure 2: -Rmatter (Eq. (5)) versus the baseline for several narrow band beams. The CP 
phase S is taken to be zero and sin'^{29i3) — 0.05. 
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Figure 3: The electron event number versus the CP phase with and without the 
matter effect. sin^(2^^i3) is assumed to be 0.05 except for the dotted curve which is for 
sin^(2^^i3) = 0.03 to show the effect of varying 6*13. Representative total errors are also 
shown. The MSD sign is assumed to be I. 
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Figure 4: Same as Fig. 3, but for different MSD signs with matter effect. 
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Figure 5: The electron event rate versus the CP phase for different narrow band 
The MSD sign is assumed to be I. 
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Peak Energy (GeV 



Figure 6: Fcp{S) (Eq. (8)) versus the peak energy of the narrow band beams. The MSD 
sign is assumed to be I. With the approximate relation, Fcp{180° + S) = —Fcp{S), the 
curves for 6 = 180°, 210°, 240°, 270°, 300°, and 330° can be inferred. 
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Figure 7: The electron event number versus the CP phase for different sin^(2^i3) values. 
The MSD sign is assumed to be I. Total errors at some points are also shown. 
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Figure 8: The CP phase versus sin^(2^i3) for a given electron event number Ne- The 
sohd (dashed) curve is for A^^e = 60 (40) at L=300 km (2100 km) with a narrow band 
beam E^peak = 0.7 GeV (4 GeV). The MSD sign is assumed to be I. 
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Figure 9: Electron event number at L=2100 km versus L=300 km for different sin^(20i3) 
values. The CP phase 6 increases from 0° (solid bullets) to 180° (crosses) then to 360° according 
to the direction indicated by the arrows. The MSD is assumed to have the sign I. In the lower 
digram for -Epeak(300) = 0.7 GeV and £'pcak(2100) = 6.3 GeV, the ellipses collapse into line 
segments. The typical total errors are also shown. 
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Figure 10: Similar to Fig. 9 for different MSD signs with fixed sin^(2^i3) = 0.05. Tlie 
results without matter effects are also plotted. 
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Figure 11: The relative error ANg/Ng versus the detector size for a 4 GeV narrow band 
beam. 
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